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ABSTRACT 

The effect of plioton-beam-induced turbulence on propagation of radio emission in 
a pulsar magnetosphere is discussed. Beamed radio emission witfi a high brightness 
temperature can generate low-frequency plasma waves in the pulsar magnetosphere 
and these waves scatter the radio beam. We consider this effect on propagation of 
radio emission both in the open field line region and in the closed field line region. 
The former is applicable to most cases of pulsar radio emission where the propagation 
is confined to the polar region; it is shown that the induced process is not effective 
for radio emission of moderately high brightness temperature but can have a severe 
effect on giant pulses. For giant pulses not to be affected by this process, they must be 
emitted very close to the light cylinder. We show that the induced process is efficient in 
the closed field line region, inhibiting propagation of the radio emission in this region. 

Key words: Plasmas-radiation mechanisms: nonthermal-pulsars: general 



1 INTRODUCTION 

It is generally thought that pulsar radio emission is pro- 
duced in the open field line region well inside the light cylin- 
der (LC) (the radius at which the corotation speed equals 
c). Radio emission generated inside the LC must propa- 
gate over a substantial part of the magnetosphere before 
escaping to the interstellar medium. The intense radio emis- 
sion is subject to various interactions with the interven- 
ing plasmas within the pulsar magnetosphere. Two types 
of interactions have been considered: cyclotron absorption 
in which electrons (positrons) are in cyclotron resonance 
with the wave (Blandford & Scharlemann 1976; Fussell, 
Luo, & Melrose 2001), and induced Compton scattering 
in which t he radio wave is scattered to a higher or lowe r 
frequency llBlandford fc Scharlem'annlll976l:IPetrovdl2004bl) . 
Both processes can affect or even disrupt the propagation 
of the beam. So far, no clear observational evidence for 
the absorption feature associated with the cyclotron reso- 
nance has been identified. However, the conventional po- 
lar cap theory, which assumes that the radio emission is 
generated in an outflowing electron-positron pair plasma, 
does predict a cyclotron resonance region where absorption 
should occur. Induced Compton scattering can be impor- 
tant in principl e for any coherent em i ssion with a sufficiently 
high intensitv JWilson fc Reejll97d ISinceU fc Kr'oiinil992t 



ICoppi. Blandford. fc Reed il993t iHardv fc Melrosd Il99, 
While cyclotron absorption is important only in the region 
near the LC, induced scattering can occur throughout the 



magnetosphere. Apart from these two processes, there is a 
third process that involves induced wave-wave interactions 
and that can significantly affect propagation of the radio 
beam. A highly collimated beam of radio waves can generate 
low frequency waves through three-wave interactions in the 
pulsar magnetospheric plasma in a way similar to wave gen- 
eration in a particle beam instability. These low frequency 
waves interact with the radio waves, leading to a decrease in 
its intensity and an increase in its angular spread. In some 
cases, such effects can completely disrupt the radio beam: 
The medium becomes opalescent to the radio waves. 

In this paper, we discuss the effect of induced 
three-wave interactions on propagation of radio emis- 
sion in a pulsar magnetosphere. Induced three- wave 
processes have been considered for pulsar winds 

iLuo fc Melrosd Il994) . pulsar eclipse b y the stellar 

wind of its companion ([Thompson et al.l 1 19941 : iMelrosj 
119941: iLuo fc Melrosdll995[r and radio emission from active 
galactic nuclei iKrishem '198^; G anga dhara fc KrishanI 
1993'; Lev jnson fc Blandford 1995). Apart from a dis- 
cussion of Raman scattering (a particular case of 
three-wave interactions, cf. Sec. 2.5) in a pulsar mag - 
netosphere JCaneadhara fc KrishanI Il993l : lLvutikovlll99ai . 
we are aware of no detailed study of the three-wave 
effect which explicitly takes account of the beaming and 
broadband nature of the radio emission. To our knowl- 
edge, there is no detailed study of propagation of radio 
emission in the closed field line region. Here we consider 
the three wave effect both in the open field line region 
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(OFLR) and in the closed field line region (CFLR) . In most 
models, propagation of radio emission is confined to the 
OFLR where the relativistic electron positron pair plasma 
streams along magnetic field lines at a relativistic velocity. 
Because electrons and positrons contribute to three-wave 
interactions equally but with opposite sign, three-wave 
interactions are significant weakened in a quasineutral pair 
plasma. However, such processes can still be important for 
radio emission with very brightness temperature provided 
that there is a modest excess of electrons or positrons and 
that the bulk Lorentz factor is modest. Since for young 
pulsars, the conventional polar cap model generally predicts 
a plasma density too high for propagation of radio waves, 
we assumes that the plasma is highly inhomogeneous across 
the polar cap and that radio emission is generated at a 
frequency near the plasma frequency (in the plasma's rest 
frame) in the underdense region where the plasma density 
is relatively low (cf. Sec. 3.1). This assumption seems to be 
required for consistency i f the radio emission is produced 
due to a beam instability jMelrose fc Gedalinl ll999'l. 

Propagation of radio emission in the CFLR is relevant 
at least for three specific cases: (1) the recently discovered 
double pulsar system PSR J0737-3039B in which the radio 
emission from one pulsar interacts with the plasma in the 
CFLR of the other, (2) radio emission from outer magne- 
tospheric regions where the emission of the trailing compo- 
nent may propagate to the CFLR as a result of the aberra- 
tion effect, and (3) backw ard emission predicted by the os- 
cillat ory polar gap model IILevinson. Melrose. Judge fc Luol 
|2QQ^ and by recent models for the pec uliar 'notch-like' fea- 
ture of pulse profiles JPvks et al J'2005). For case (1), eclipse 
of the radio emission from one pulsar by the magnetosphere 
of the other has been observed. A possible interpretation of 
the eclipse is due to processes in the CLFR such as cyclotron 
or synchrotron absorption or induced scattering. For case (2) 
the trailing component may have an absorption-like feature 
or may be completely destroyed as the result of strong in- 
duced processes. For the third case, induced processes can 
severely constrain visibility of the backward emission that 
propagates through the CFLR. 

In Sec. 2, we outline a general formalism of induced 
three-wave interactions in the pulsar magnetosphere. We 
discuss the effect on propagation of radio emission in the 
OFLR in Sec. 3 and in CFLR in Sec. 4. 



2 INDUCED THREE- WAVE INTERACTIONS 

We consider three- wave processes due to passage of beamed 
radio emission (referred to as a photon beam or radio beam). 
The processes involve two high frequency waves (unscat- 
tered and scattered radio waves) and a low frequency wave 
generated by the radio beam. The three waves satisfy the 
resonance conditions w w' -I- uj" and fe ^ fe' -I- fc". The 
condition u ~ oj' ^ ui" means that three wave scattering is 
approximately 'elastic'. 



2.1 Three- wave probability 

We model pulsar radio emission as a highly coUimated 
photon beam, described by the photon occupation num- 
ber N{k), where uj is the wave frequency and k is the 



wave vector. The photon occupation number is a Lorentz 
invariant, related to the brightness temperature Tb through 
hu>N{k) — ksTb, where fes is the Boltzmann constant. Let 
NL{k") be the occupation number of the low-frequency wave 
with frequency u)" and wave vector fc". Evolution of N{k), 
N{k') and NL{k") are determined by three kinetic equations 
which are completel y determined by the three- wave proba- 
bility jMelrosell986ll : W{k, k'k"){2TvfS{k-k' ~k")5{Lu{k)- 
Lu'{k') — Lo" (k")), where the delta functions are the three- 
wave resonant conditions and 



W{k,k',k") = TiRl 



(1) 



In Q, which is simply referred to as the three- wave prob- 
ability, Uija is the quadratic response tensor, Rl is the 
ratio of the electric to total energy of the low frequency 
wave, and the polarizations of three waves are e, e' and e" . 
In deriving the three-wave probability, we adopt the cold 
pl asma form for the qua dratic response tensor, e.g. Eq (Al) 
in lLuo fc Melrosel ll2006f) . and take the strong magnetic field 
limit. One should note that the strong field limit may not 
be appropriate for the region near the LC. Nonetheless, it 
allows us to simplify the quadratic response tensor greatly. 
In this limit, one finds that the dominant component is the 
projection to the magnetic field, denoted by a«, given by 



rjEoeu} , 
ail « n\\ 







(2) 



where ui u' ^ lo" , uop is the plasma frequency, ny = 
k\\c/uj, n'l = k'^^c/u', nj' = fcj'c/a;", rj is the ratio of the 
difference to the sum of the positron and electron number 
density, with r] = —1 corresponding to an electron gas. We 
use the subscript || to denote the parallel (to the magnetic 
field) component. The final expression for the three- wave 
probability is 



vr 2 2 ft'^ /t^' 

W ^ —rj reC ; 

2 nieC 



(3) 



where = /{i-KSomeC'^) ~ 2.8 x 10~^^ m is the classical 
electron radius and = {cup/uj")^^^ {n^^ + n'^\ — n'0e\\e'^^e'^' . 
The expression ^ is closely anal ogous to a sim ilar expres- 
sion for an unmagnetized plasma iMelroselll994r) where the 
low frequency waves are the Langmuir waves {u;" = Up). 
For a quasi-neutral, cold electron-positron plasma, one has 
?7 = 0; one then needs to retain higher order terms in an ex- 
pansion in {uj,uj")/fle ^ 1 to obtain a non zero three-wave 
probability, which is much smaller than iLuo fc Melros j 
1997). 

In the following, we consider two low frequency modes 
that may be relevant to induced three- wave interactions: the 
Langmuir mode that propagates along the magnetic field, 
and the Alfven mode (also called the low frequency LO 
mode) which can propagate at any direction. Note that the 
three- wave process involving a Langmuir wave in the (coher- 
ent) fixed-phase formalism is also called stimulated Raman 
scattering (cf. Sec. 2.5). Landau damping imposes a lower 
limit to the wavelength of low frequency waves. To deter- 
mine the Landau damping one need consider separately an 
intrinsically relativistic plasma such as in the OFLR and a 
nonrelativistic plasma such as in the CFLR. Landau damp- 
ing in an intrinsically relativistic plasma occurs preferen- 
tially for n" ^ 1 since the plasma is relativistic even in the 
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plasma rest frame and the particle's Cerenkov resonance oc- 
curs for the wave phase speed near the speed of light. In 
the nonrelativistic plasma, Landau damping is effective for 
electrostatic waves with n" ~ w^^^ = 0.3c/AotJp (where 
Ad is the Debye length) iBekefil |l966l , leading to a limit on 
the maximum n'^^y^ ^ 1. Since l^pj^ oc n"'^, three-wave in- 
teractions involves low frequency electrostatic waves should 
be important in the CFLR where the plasma is expected 
to be nonrelativistic (cf. Sec. 4.1). In the OFLR the plasma 
is intrisically relativistic, Landau damping limits the refrac- 
tive index of the Langmuir waves to about 1, and both the 
Alfven and electrostatic modes need to be considered. 



2.2 Scattering rate 

We consider the production rate of scattered photons as a 
measure of the effect of the three-wave process on the pho- 
ton beam. This rate, which we also call the scattering rate, 
depends on the level of the low frequency waves generated 
by the photon beam. We consider the case where the low 
frequency waves are saturated as the result of change in the 
occupation number of the photon beam. The saturation con- 
dition dNL/dt — leads to a level of low frequency waves 
expressed in terms of the photon occupation number. Thus 
the scattering rate can be written as 



r = 



dlniV(fc') 
di 



Af7P WN(k) 



(4) 



(27r)2 c ' 

where Afl is the angular spread of the photon beam. Our 
specific model for an axially symmetric beam is 



iV(fc) = N{k)b{e), 6(6) = exp - 



2Aeg 



(5) 



where Af2 = 27rAOo, AOo <g 1 is the beam width and the 
polar angles (O,"!?) are defined with respect to the beam 
axis. Similarly, we use (0",$") for the polar angles of k" . 
The polar angles of k and k" relative to the magnetic field 
{z axis) are given respectively by (0,4)) and (6" ,(f)"). Since 
the angular spread of the beam is small, hereafter we do not 
distinguish between the beam axis and fe. Substituting ^ 
for Q yields 



(6) 



Stvc V rUeC^ 

where we express the incoming photon occupation number 
in terms of the brightness temperature: TtwNik) — ksTi,, 
and depends on the specific low frequency mode involved 
and is derived in Sec. 2.3 and 2.4. 

From @ one may define the opacity for the photon 
beam propagating over a distance As: 

^As 



(7) 



In ||7|l, AfiTi, can be estimated from the observed fiux den- 
sity SZ^"": 



Ann 



1/2 \AsJ 



V 105 my V 



IGHz 



Dl 



1 kpc 



lOmJy 



(8) 



where = uj/2n is the frequency in Hz, As is the distance 
from the emitting region to the region concerned and Dl is 
the pulsar distance. If one assumes the size of the emission 
region is 51, one has Afi ~ {31 /As) . The brightness tem- 
perature is estimated to be TJ, '-^ 5 x 10^* K for SI = 10^ m. 

2.3 Back scattering 

Consider first the case where the low frequency waves are 
Langmuir waves propagating along the magnetic field at a 
low phase speed uj"/k" -C c (or n" S> 1). Three- wave in- 
teractions favor the back scattering regime where the in- 
coming and scattered photons propagate in opposite direc- 
tions. From the three-wave resonance conditions one has 
fc|| = -I- k" (parallel to the magnetic field) and k± = fe'^ 
(perpendicular to the magnetic field), which gives rise to 
cosO ~ cos^' -f- n" [ijjp/ uj), where n — kc/uj ^ n = k'c/u' ~ 
1 and uj" — uip. Because of fc « k' , one must have ~ vr — 
corresponding to the incoming and scattered photons prop- 
agating in opposite directions. Hence for Langmuir waves, 
one derives 



//2 . 4 , 

n sm ( 



(9) 



provided that 2 cos 6'(aj/tjp) < nj^ax; where nj'jax is de- 
termined by the Landau damping. The angular range for 
back scattering to occur is cos 61 ^ mm{l, k'^^^/2k} with 
^max = 'T-max'^p/c. When the photon beam is parallel to the 
magnetic field, one has |^pp ~ (tj/ajp)^(Af2/7r)^, which im- 
plies that only photons off the beam axis with 6* 7^ can in- 
teract with the Langmuir waves. In a nonrelativistic plasma, 
one has nj^ax ^ 1 and hence the back scattering can be very 
efficient. 



2.4 Diffusion of the photon beam 

In the regime k" <^ k k' , the three- wave interactions re- 
semble particle- wave interactions, with high frequency waves 
(the radio emission) acting li ke 'parti cles' that can emit or 
absorb low frequency waves jMelrosei ri994: Lu o fc Melros^ 
Il995l [l997l) . This regime is also called the small angle scat- 
tering because the angle between the incoming and scattered 
photons is small, given by 6'^.^,/ ~ k" /k <^ 1. The three- wave 
resonance conditions reduce to the Cerenkov resonance con- 
dition between 'particles' (the high frequency wave) with a 
velocity Vg = duj{k)/dk and a wave with [uj" ,k"). Specifi- 
cally, one may define a Cerenkov angle xo = arccos((jj"/fc"c), 
corresponding to a conical surface with respect to Vg, which 
defines the Cerenkov resonance condition. The condition can 
only be satisfied for a subluminal wave n" = k" c/uj" > 1, 
i.e. its phase speed must be less than c. We are interested in 
negative absorption (or growth) of the low frequency waves, 
which occurs for Q" ^ Xo ~ AOo, at a rate determined b y 
the three- wave probability Q or ||3J| iLuo fc Melrose|[l997ll . 

The instability of low frequency waves leads to diffusion 
of photons in the fe-space, determined by the diffusion coef- 
ficient Dij oc k"^ k" kjW Nl /c, where the three-wave proba- 
bility W is given by Q . The growth of low frequency waves 
is entirely due to the anisotropic angular beaming of the 
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photons. A decrease in the beaming, over the perpendicular 
diffusion time t± = k^/D = AQk^/D, where D = Da, 
is especially important as it can ultimately suppress the 
growth. Assuming that the low frequency waves are sat- 
urated due to perpendicular diffusion, one can derive the 
total diffusion time td = k^/D = l/F, which describes the 
effect of the instabihty on th e beam. One finds F is the same 
as that given by Q or © JLuo fc Melrosdll995l) and the 
diffusion effect can be estimated by evaluating the opacity 
0. For Langmuir waves, one finds 

^ (A cose -^^ sine -iV 

X [^1 + I^-J^ cot e^ n"'sin^e, (10) 

where 1 < n" ^ Lojujy. Similarly for Alfven waves one has 

(2cos6»- 1)' sin'' Stance", (11) 

where the parallel component of the polarization is 
le ^l ^ (tj"/a;„)^n''sine" < 1, d" < 1 and nj,' ^ 
1 JArons fc BarnarJigsdiMelrose fc Gedalinlll999ll . 



have an inverted energy population, and this has b een sug- 
gested as a mechanism for generating giant pulses dPetroval 
|2004a*). The induced three-wave interaction discussed here 
is essentially the same process described as 'stimulated Ra- 
man scattering' by Gangadhara and Krishan (1993) and 
as 'induced Raman scattering' by Lyutikov (1998) except 
that in our case the low frequency wave is not restricted 
to plasma waves. The important distinction, as noted by 
Lyutikov (1998), is that in the induced three- wave interac- 
tion, the Raman frequency, ujr, is replaced by the frequency 
of the low frequency waves, which is the plasma frequency, 
LOp, in cases of interest. However, this analogy is only kine- 
matic, and the dynamics of the two processes are quite dif- 
ferent; for example, they are described by different kinetic 
equations. In the induced three- wave interaction the plasma 
may be regarded as passive, determining the wave dispersion 
and the nonlinearity that allows the three-wave coupling. 
The kinetic equation involves on the wave distributions: the 
low frequency waves grow and result in the scattering of 
the high-frequency waves. In contrast, in induced Compton 
scattering, the plasma particles actively scatter the high- 
frequency waves, the process depends explicitly on the dis- 
tribution function for the particles, and no low-frequency 
wave is involved. 



2.5 Induced Compton scattering 



The relation between the induced three- wave interaction dis- 
cussed here and induced (or stimulated) Compton scatter- 
ing requires comment. (We use 'induced' and 'stimulated' 
interchangeably here.) Both processes need to be consid- 
ered in discussing the propagation of radiation with very 
high brightness temperature, and we regard them as inde- 
pendent of each other. In our view, the difference between 
these two processes has become somewhat obscured by dis- 
cussions of stimulated Raman scattering. Formally, Raman 
scattering is associated with particles (atoms or molecules) 
that have an intrinsic transition frequency, lor say, associ- 
ated with them. Spontaneous Raman scattering of incident 
radiation with frequency cu produces scattered radiation at 
the Stokes frequency, cu — lur, and at the anti-Stokes fre- 
quency, u! + lor. Stimulated Raman scattering is analogous 
to absorption of the beat between the incident and scattered 
radiation, and occurs at a rate that is proportional to the 
product of the intensities at the incident and the scattered 
frequencies. This leads to exponential growth at the Stokes 
frequency, and an associated exponential decay at the inci- 
dent frequency, so that radiation is transferred from higher 
to lower frequency at an exponentially increasing rate. (If 
the scattering particles have an inverted energy distribution, 
negative absorption or maser action results in an exponen- 
tially increasing transfer from the incident frequency to the 
anti-Stokes frequency.) 

Induced Compton scattering is analogous to induced 
Raman scattering, except that there is a continuum of fre- 
quencies for the scattered radiation. For isotropic radiation 
with a peak in the intensity at some frequency, induced 
Compton scattering tends to move this peak to lower fre- 
quencies; if there are two distinct beams of radiation, in- 
duced Compton scattering tends to transfer radiation from 
the beam at the higher to the beam at the lower frequency. 
The direction of transfer is reversed if the scattering particles 



3 PROPAGATION IN THE OFLR 

The formalism derived in the preceding section can be ap- 
plied to the case of propagation of radio emission in the 
OFLR where the plasma has a bulk velocity v — f3c and a 
Lorentz factor 7 = 1/(1 — fd'^Y''^. A convenient approach is 
to treat the formalism in Sec. 2 in the plasma rest frame and 
then derive the relevant quantities in the pulsar or observer's 
frame by a Lorentz transformation. Note that some caution 
is needed in defining the bulk rest frame. Since the plasma 
streams along curved magnetic field lines, there is no single 
bulk frame. However, one can still construct a consecutive 
series of local bulk frames along the curved field lines, an ap- 
proach similar t o that used in ray tra cing in a relativistically 
flowing plasma llFussell fc Ludl2004f) . We assume that the 
photon beam is axially symmetric in the plasma rest frame. 
The photon beam shape varies from frame to frame. Such 
complication is not considered here, though it should be in- 
cluded in a full radiative transfer calculation. Specifically, 
the relevant propagation angle is transformed according to 



h sin 9 • 



n sin 9 



7(1 — n(3 cos t 



n cos o — 



n cos 9 — P 
nl3 cos 6 



(12) 



All tilde quantities correspond to the rest frame. We ig- 
nore the pulsar's rotation in the Lorentz transformation. 
The brightness temperature, frequency, differential solid an- 
gle and plasma frequency are transformed according to 
fb = 7(1 - nl3cos9)Tb, = 7(1 - n/3cos9)uj, AQ = 
Af2/7^(1 — nf3 cos 6)^ , and Up = uipf-y^^^, where the refrac- 
tive index is n ~ 1. Since the time derivative in Q is inter- 
preted as the total time derivative, i.e. d/dt = d /dt + cn ■ V 
with n = k/k, one has f = r /7(l — 0cos6) I^ (j)mr^nin3 
119731: iHardv fc Melrosd Il995l: ISazonov fc SunvaevI Il99'^ 
Upon making the Lorentz transformation one obtains the 
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scattering rate in the observer's frame: 

^ retype / A^fcaTh X (^pe 

87r(2Af7)i/2c V mec2 J 1 + 6* V ' 

where M is the pair multiphcity. In (1131 . we use r] — 1/2M 
and rewrite the plasma frequency as u>p = {2My''^ujpE, 
where Upe is the plasma frequency corresponding to the 
number density of the excess electrons or positrons. 



approximation provided that the size the scattering region 
(r,/7) is much larger than the length scale (corresponding 
to the polar gap length) of plasma inhomogeneities along the 
propogation path. These two assumptions allow us to esti- 
mate (I13II without going into much details on a particular 
polar cap model. We consider the following three propaga- 
tion regimes separately: two quasi-parallel regimes, includ- 
ing the forward propagation 9 <^ 1 and backward propaga- 
tion S ~ TT, and quasi-perpendicular regime 6 ~ n/2. 



3.1 Plasma density 

The scattering rate is proportional to the local net 
charge density: a higher net charge density implies a 
faster scattering rate. This is due to the electrons and 
positrons contributing with opposite sign to the three- 
wave coupling coefficient. For the conventional polar cap 
model the number density of the excess charge is the 
Goldreich- Julian (GJ) density, given by ngj « 7 x 
10^m-3(P/0.1s)-i(S/10*T)(7?o/r)3 dColdreich fc JuhanI 
Il969h . where Ro — 10* m is the star's radius. Seemingly 
insurmountable difficulties arise if one assumes the most 
widely discussed emission mechanism (a beam instabil- 
ity) in a plasma with the GJ density times a multiplic- 
ity M ^ 1. The lowest frequency consistent with theory 
is (jJc ^ 2{'y/{i)Y/'^ujp jMelrose fc Gedalinlll999l) . where 
{7') is the Lorentz f actor's mean spread in the plasma rest 
frame, about 1-10 iArendt fc Eilekl l2002h . In the follow- 
ing we assume {7') = 4 implying loc = ^^^'^LOp. For ex- 
ample, the plasma frequency predicted for the Crab pulsar 
is tjp « (2 X 10^^s-^)(i?o/?■)^''^(M/100)^''^ corresponding 
to ~ l^^'^iOp f« (2 X 10^^ s-^){Ro/rf''^{M/lQQf''^ for 
7 = 100, which is too high compared with the observed 
frequency. (This plasma density problem was discussed by 
[KuuzI. Lcsch, Jcssrieii& Hocr^^^^ 1.1998.) .) One practi- 

cal solution to this problem is that the plasma is highly in- 
homogeneous forming underdense regions where radio emis- 
sion is produced and propagates. The low density region 
where the density is actually much lower than 2MnGj 
can form as a result of nonuniform extraction of pri- 
mary particles from the polar cap or nonun iform pair cas- 
cades iLevinson. Melrose. Judge fc Luol2005^ . In the former 
case, the primary particle density is nonuniform across the 
polar cap forming localized outflowing fluxes along field lines 
separated by low density regions where the net charge den- 
sity is much less than ncj (though on average the mean 
flux density from the polar cap is still nojc). Both the mul- 
tiplicity M and the net charge density can be low in the 
underdense region. 

In the following we consider the case where the ra- 
dio wave is emitted in the low density regions at a fre- 
quency near the local plasma frequency. Specifically we as- 
sume (1) the radio frequency is ct; ~ ujpt,^^^^ , where ujp* ^ 
{e^MnGj/£ome)^^'^{Ro/rf)^^^ is the local plasma frequency 
at the emission radius r* and M is the mean multiplicity, 
and (2) AI ^ M is a constant along the propagation path. 
Assumption (1) is applicable if the c oherent radio emis- 
sion i s due to wave-particle interactions iMelrose fc GedalinI 
^23). Although assumption (2) may not be realistic, par- 
ticularly in the oscillatory case where both the multiplic- 
ity M and the net charge density ca n be nonstationary 
iLevinson. Melrose. Judge fc Luoll2005l) . it is a reasonable 



3.2 Quasi-parallel propagation 

In the forward propagation regime 6 <ti 1, fc is nearly parallel 
to the magnetic field. Assuming the low frequency plasma 
wave propagates along the magnetic field, one has 6" = 
and h" ~ 1. We assume that the photon beam is in the 
LO mode with the polarization in the B-k plane, because 
only this high frequency mode can be produced through 
wave-particle interactions jMelrose fc Gedalinlll999l) . If one 
assumes that the radio emission is generated near the fre- 
quency Lu ~ cjp»7^''^ at the emission radius r«, the plasma 
frequency in the rest frame is comparable with the photon 
frequency. Thus the only relevant low frequency waves are 
Alfven waves. 

For the Alfven mode, one can show that three-wave 
interactions due to the projection of the polarizations of 
all three waves along the magnetic field are small. From 
mi l, which applies in the rest frame, one obtains |^p|^ ^ 
{Co" /Qpfe'' < 1, where e\\ ~ e\ ~ 9. For Cj" /up = 0.2 and 

9 = 0.1, one has ||pp ^ 10"^ 

In most part of its propagation path the radio emission 
propagates backward in the rest frame, as indicated in fig- 
ure0 and in this regime the low frequency approximation is 
justified. From lO, one has cos^ ^ (1 - 6l^7^)/(l -I- 6*^7^) ^ 
— 1 for 9 3> 1/7, i.e. the emission propagates backward 
(9 ~ tt) in the rest frame. In the rest frame the fre- 
quency is given by a) ~ 9'^^uj/2. Since lj ~ iOpt^^^'^ , one 
has uj/Cjp ^ 0^7^ S> 1 or S> oJp, where tip ^ ujp, for 
As/r, = {r — r,)/r, < 1. Thus for backward propaga- 
tion, the Langmuir waves can be treated as low frequency 
waves (ti ^ Ljp) and back scattering (cf Sec. 2.3) is allowed 
kinematically but is inhibited by the Landau damping. For 
backward propagation, the LO mode is elliptically polar- 
ized, with only a small component along the magnetic field 
e|| I ~ |e|| I « 2/9'Y ^ 1. From the three- wave resonance con- 
ditions one finds n" ~ 2Gj jihp ~ 61^7^ 3> 1. However, Landau 
damping limits the refractive index to n" ~ 1. 

For backward propagation, the three-wave interaction 
can proceed in the regime fc" ^ fc ~ fc' where both Lang- 
muir and Alfven waves are relevant. From llUl one obtains 
2|IpP/(1 + ^'^7^) ~ 2/6''^7'' < 1 for the Langmuir waves with 
n" « n" ~ 1 and from one has 2l4P/(l + ^^7^) ~ 
2{Cj" /Cbpf{9"f /9'^-i^ for the Alfven mode. Thus, Alfven 
waves are less efficient than Langmuir waves. 

3.3 Quasi-perpendicular regime 

The transition between forward and backward propagation, 
9 = 7r/2, in the rest frame corresponds to 5 ~ I/7 in the 
pulsar frame. A photon beam that is emitted initially along 
the magnetic field rapidly acquires an angle as the magnetic 
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Figure 1. Propagation angle in tlie rest frame vs propagation 
distance As/r«. The propagation angle 6, wliicli is assumed to 
start from zero (at As = 0), rapidly increases to tt after a distance 
As = r — r* > r» j^. A dipole magnetic field is assumed. The solid 
and dashed lines correspond to 7 = 100 and 50, respectly. 
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neutron star 




Figure 2. Propagation of the photon beam emitted at r*. The 
propagation enters the quasi-perpendicular regime at T where the 
propagation angle in the pulsar frame is S ~ I/7. The distance 
from the source (5) (at r*) to the quasi-perpendicular regime (T) 
is approximately r»/7. 

field lines bend away (cf. figure 0. The propagation enters 
the quasi-perpendicular regime [Q ~ 7r/2) after the distance 
As ~ I"*!!- Similarly to forward propagation, the plasma 
frequency is comparable with the radio frequency in the 
rest frame. So, the plasma waves cannot be treated as low 
frequency waves. The relevant low frequency waves are the 
Alfven waves. From we have ||pp {uj" j u^,)^ Q" . 
Although the process favors a higher cij", the Alfven mode 
is heavily damped near the plasma frequency. It can freely 
propagate in the range tZ) <C (ip. For example, one may 
choose B" = 0.5 and ti'Vcip = 0.2. This gives 10"^. 

To calculate 1131 1. we determine the plasma frequency 
in the propagation path from uj^e ~ iUp/{2M)^^^ — 
[a;/(2M7)^/^](r*/r)^/^ with r ~ r* in the quasi- 
perpendicular region. Note that for young pulsars such as 
the Crab pulsar the local plasma frequency determined in 
this way is much lower than the one derived from the den- 



sity MnGj{Ro/r,)^ (cf. Sec 3.1). The opacity r = FAs/c = 
Fr«/7C is estimated as 




where JHJ with As — r,/7 is used. It should be noted that if 
the emission is from a region spread over r*/7, the opacity 
estimated from 1141 is reduced. If one assumes the nominal 
value I^pI^ — 10^^ for typical pulsar parameters, Eq. 1141 
shows that the induced diffusion is not important for radio 
emission with a moderately high brightness temperature. 
However, it can be effective for young pulsars or pulsars 
that emit giant pulses. 



3.4 Giant pulses 

For giant pulses the opacity 1141 can well exceed 
unity and as a result, it imposes a severe con- 
straint on their propagation. Giant pulses have been 
observed from a few pulsars wi th the two be s t ex- 
amplcs being the Crab pulsar ("Lundg ren et alJ Il995l : 
[llankins, Kern, Wcathcrall. & Eilek 2003 ), a fast rotating 
young pulsar and PSR 1937-1-21 llKinkhabwala fc Thorsetti 
l2000l) . the second fastest millisecond pulsar. The distribu- 
tion of pulse energy (the fiux density times the typical pulse 
width) is power-law, in contrast to normal radio pulses that 
have a log-normal distribution IRo mani & Johnston 200 ij). 
This implies that the emission mechanism for giant pulses 
is distinct from that for normal radio pulses. Although the 
emission mechanism for giant pulses is not well understood, 
recent study of giant pulses in relation to the corresponding 
high energy emission suggests that the emission region may 
be in or near the outer magne tospheric region where th e 
high energy emission originates llRomani fc Johnstonll20()lh . 
Here we show that the opacity constraint 1141 excludes the 
possibility of giant pulses originating from the inner mag- 
netosphere and deduce that they must come from the outer 
magnetosphere, near the LC. 

The opacity 1141 decreases with increasing emission ra- 
dius. This feature can be understood as follows. The ef- 
ficiency of induced three-wave processes depends on the 
photon energy density. For an emitting source of size 
51, the photon energy density in the quasi-perpendicular 
region (at a distance t-«/7 from the source) is propor- 
tional to (Sl/rt)^ ksTi,. Observations of nanosecond giant 
pul ses from the Crab suggest a source size o f order a me- 
ter faankins. Kern. Weatherall. fc Eilekl2003l) . This density 
decreases with increasing r* , leading to a weakening of three- 
wave interactions. For giant pulses from the Crab pulsar, one 
has an average flux density 5'°'^= ~ 10^ Jy. We have r > 1 
unless a large emission radius r« is assumed. For giant pulses 
not to be strongly affected by the induced process, corre- 
sponding to the condition r < 1, the emission radius must 
satisfy the condition 
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where \^p\'^ = 10~^. Although the pair multiplicity is not 
well constrained, HSU still gives a strong constraint on the 
emission radius for a reasonable parameter range. A re- 
cent numerical simulation suggests that the conversion ef- 
ficienc y of the primary parti cle energy to pairs is well below 
unity iArendt fc Eilekll2002^ . corresponding to M7 < 76/2, 
where 75 is the primary particle's Lorentz factor. For young 
pulsars, resonant inverse Compton scattering is the domi- 
nant emission process for pair-producing photons and it is 
reasonable to assume 76 ~ 10^ This gives an upper limit 
M7 <S 5 X 10*. Applying 1151 to the Crab pulsar, assuming 
a distance Dl = 2 kpc, M < M = 10^ and 7 — 10^ one finds 
> 5x10^^ m, about 0.3 of the LC radius _Rlc ~ 1.6x10*^ m. 
One concludes that for giant pulses to propagate without be- 
ing dispersed by induced three-wave interactions, they must 
be emitted very close to the LC radius. For the millisecond 
pulsar PSR 1937-1-21, one expects M to be small and 7 to be 
large as pair production in the millisecond pulsar's magnetic 
field is much less efficient. Assuming that the primary par- 
ticle's energy is limited by curvature radiation correspond- 
ing to 76 ~ 10^, one has the upper bound M7 < 5 x lO''. 
For example, for Dl ~ 3.6 kpc, M = 10, 7 = 10^ and 
5°*"" = 300 Jy, one has r* > 3 x 10'' m, about half the LC 
radius. 



4 PROPAGATION IN THE CFLR 

The formalism discussed in Sec. 2 can be directly applied to 
the CFLR where the plasma consists of single-sign charged 
particles with the GJ density ncj and is stationary with a 
zero bulk velocity. All the relevant quantities discussed in 
Sec. 2 are now referred to the pulsar frame. 

4.1 Plasmas in the CFLR 

In contrast to the polar region, the plasma in the CFLR is 
plausibly stationary in the corotating frame. Since particles 
are in the ground Landau state they do not have a per- 
pendicular momentum, so there is no magnetic mirror effect 
that can trap a dense plasma. This feature is very different 
from the usual planetary magnetosphere in which a plasma 
is trapped due to this effect. Thus, in the absence of continu- 
ous injection of pairs the plasma density in this region is con- 
trolled by the pu lsar's rotation in such way that the charge 
density is encj iColdreich fc Julianlll969fl . The plasma in 
this region consists of particles with a single sign of charge 
(either electrons or positron or possibly protons), which can 
be treated as an electron gas in one dimension (along the 
field lines) . Assuming the magnetic field as a function of ra- 
dius r: B on {Ro/rY , where (5 = 3 for a magnetic dipole, one 
may express the density as Ue = ncj {Ro/rY. Using and 
J^, one obtains 

~ . f. ( Anr6 \ / u y ( P 

^ ^ V3 X 1022 k; Vl GHz/ V50msy 

nio^J [^—) (0:2^)1^-1 ' 

where is given by (|5J for back scattering and 1101 for the 
quasilinear diffusion approximation. Notice the difference in 
the radial distance dependence in 1161 compared to 11411 . 
This difference is due to the GJ density varying with radial 



distance proportional to the magnetic field, where we ignore 
a cosine factor. The CFLR region is separated into a positive 
charge region and a negative charge region, with the bound- 
ary between them called the null surface. One may define a 
radius r^r where r = 1. The opaque region corresponds to 
r ^ Vcr, where induced three-wave processes are important. 
This critical radius is derived as 




We are interested in the region sufficently far from the star's 
surface such that the condition uj S> (j->pe is satisfied. We 
discuss only the case where the low frequency waves are 
electrostatic. 



4.2 Eclipse in the binary pulsar J0737-3039B 

We first apply H16^ to the double pulsar system J0737- 
3039B, in which radio emission from one pulsar is period- 
ically eclipsed by the other. It is of interest to consider if 
induced scattering can provide an effective mechanism for 
such eclipse. The system consists of pulsar A with the pe- 
riod 22.7 ms and magnetic field 6.3 x 10^ T, and pulsar B 
with the period 2.8 s and magnetic field 1.6 x 10** T. The two 
pulsars are in a nearly circular orbit (with the eccentricity 
0.088), with the orbital plane seen nearly edge-o n, at an in- 
clinat ion angle ~ 87° to the line of the sight (iLvne et alJ 
l2004h . Observations indicate that the nearest impact dis- 
tance of the photon beam from pulsar A to pulsar B is about 
O.IRlc — 0-2Rlc (of pulsar B). The plasma frequency at 
r = O.IRlc is about Upe ~ 2 x 10^ s~^. Cyclotron absorp- 
tion is p robably unimportant because the plasma density is 
too low jLvutikov fc Thompsonlboo'sll . 

Consider interaction of a photon beam from pulsar A 
with the magnetospheric plasma of pulsar B. The brightness 
temperature of the radio emission from pulsar A is estimated 
to be Tt - (4.8 x 10^^ K)/An at 1.36 GHz at r = O.IRlc « 
1.4 x 10^ m. To estimate the opacity, one may take r ^ 
As = O.IRlc- Landau damping limits the Langmuir waves 
ton" < n'^^^ « 2xlO^(Te/10^ Ky^^^{n^/2.5xl0'' m^'-^Y^^ , 
where Te is the plasma temperature. For P = 2.8 s, and 
n" = 2 X 10^, one estimates r ^ 10'"'. The possibility of 
injection of relativistic pairs in the magnetosphere of pul- 
sar B was recently discusse d by Rafikov fc Goldrcich (200£f), 
also lLvutikov fc Thompson! ll2005l) . In Rafikov fc Goldreich's 
model, trapping of a dense pair plasma is possible due to 
particles' nonzero pitch angles that are acquired as a result 
of synchrotron absorption of radio emission from pulsar A. 
Since the opacity is proportional to the net charge density 
and is unchanged by the addition of pairs. Thus, we con- 
clude that the induced three-wave process is not important 
for this binary system. 

4.3 Asymmetry in double-peak profiles 

Induced three-wave processes may be relevant in the outer 
magnetosphere for radiation originating near the last open 
field lines. Waves can be swept into the CFLR and this ef- 
fect is especially important for a class of pulsars with double 
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peaks with a wide separation of nearly 180°. In the single 
pole model, the main (the leading main) and inter (the weak 
trailing) pulses originate from emission from the opposite 
sides of a wide cone at rather high altitudes. Two compo- 
nents are often asymmetric in intensity with one (normally 
the leading component) being stronger than the other. One 
plausible interpretation of such asymmetry is that waves 
from the trailing component are partially absorbed due to 
cyclotron absorption (.Fussell fc Luo 2004') . Since the trailing 
component can be swept into the CFLR, induced three- wave 
process is a possible alternative mechanism for its suppres- 
sion. 

4-3.1 Radio beam crossing the CFLR 

To derive the condition for the radio emission to propagate 
across the CFLR we assume that the photon trajectory is 
straight in the observer's frame. In the co-rotation frame, 
the photon trajectory appears as a spiral curve. We con- 
sider an orthogonal rotator with the line of sight crossing 
the magnetic axis in the equatorial plane and emission from 
one pole only. In the corotating frame, one assumes a static 
dipole with field lines described by r = C,^Rlc s\v? ^P, where 
C, — sin 'I'd/ sink's, 'I'd is the half opening angle of the po- 
lar cap, and 'ts is the colatitude of the field line foot point 
on the star's surface. The last open field lines correspond 
to = 1. Assume that waves are emitted at the radial dis- 
tance r, on the last open field field lines. After propagating 
a distance s the radial distance is approximately r, + s. If 
the photon beam propagates along a straight line given by 
(r, -|- s,ct) in the observer's frame, where 7r/2 — cr is the 
polar angle (relative to the magnetic axis) of the propaga- 
tion direction that differs from the tangent to the field line, 
the condition for waves to propagate into the CFLR is (Ap- 
pendix): 

r,+s^('^RLCsin^{'iP + s/RLc), (18) 

where 'i' = + s/Rlc and ^ ^ 1. Note that can be 
expressed in terms of r, , s and a using Eq and lIXoll . 

This condition is shown in Figure 01 The condit ion can easily 
be satisfied for emission originating from a large emission 
radius and a field line close to the last open field line ~ 1. 




0.1 0.2 0.3 0.4 0.5 0.6 
s/Rlc 



Figure 3. The condition for a wave propagating into the closed 
field lines. The contours (from innermost) correspond to closed 
field lines of C = 0.995, C = 0.999, and C = I. The emission 
region is assumed to be on the last open field lines (f = 1). The 
dashed line represents the propagation path of a wave emitted at 
r* = 0.2Rlc- The wave enters the CLFR at s = 0.04Rlc and 
exits at s = 0.57 Rj^c- 

propagate into the CFLR. From figure this radius must 
be in the range r« < 0.13-Rlc or r* > 0.5Rlc- Since the LC 
radius is only Rlc ~ 7.2_Ro and r, > _Ro ~ O.URlc, the 
only permitted radial range is r* > 0.5Rlc- Alternatively, 
one may assume the radio emission occurs in a region away 
from the last closed field line surface and in this case it prop- 
agates only in the OFLR. However, one needs to assume a 
two-pole model to explain the double peak feature. 

For the Crab pulsar, the observed fiux density is S°^^ ~ 
950 mJy at 400 MHz and 5°''" ^ 14mJy at 1.4 GHz, cor- 
responding a spectral index a ~ 3.4. This gives A^lTt 
3 X 10^" K for i/ = 1.4 GHz at r = Rlc « 1-6 x 10^ m and 
Vcr ~ 0.5-RlcI^pI'*'"^. Since the propagation in the CFLR is 
nearly parallel, one has back scattering with |^pp ~ n"^ ^ 
1. Since both the leading and trailing pulses are observed, 
the emission radius must be in the range r« < O.ISRlc or 
r. > 0.5Rlc- 



4.3.2 Fast pulsars 

The fate of the radio emission that propagates into the 
CFLR is determined by the opacity llfcill . We consider PSR 
1937+21 and the Crab pulsar. For PSR 1937+21, we have 
the mean flux density S"*"^ « 240 mJy at 400 MHz and 16 
mJy at 1.4 GHz. The brightness temperature times AQ at 
Rlc is estimated to be AD.T1, ~ 10^*^ K and 5 x 10^^ K, 
respectively. The maximum depth from the last closed field 
lines boundary is A^niax ~ {r,+s)Aip = (r, + s)^''^/(l — r, — 
s)i/^AC ^ 2.5 X IO-^'Elc for r. + s = 0.5, AC = 0.01. For 
P = 1.5 ms one has A^max ~ 2 x lO^m. Using 11711 and as- 
suming B — 10^ T and As = 0.2Rlc, one finds r^r > Rlc- 
The opacity at r ~ Rlc is r ~ 53|CpP > 1 at = 1.4 GHz 
and it is even larger at a lower frequency. This implies that 
propagation in the CFLR is inhibited. Since both the lead- 
ing and trailing components are visible, the emission must 
be produced at a radius where its trailing part does not 



4.4 Constraint on backward emission 

Induced three-wave processes impose a severe constraint on 
any model involving backward emission in the pulsar frame. 
Backward emission may b e relevant for oscillatory gap mod- 
els su ch as proposed by iLevinson. Melrose. Judge fc Luol 
(|200^. In this model, pair cascades occur in an oscillatory 
manner sending secondary particles both forward and back- 
ward. This would lead to two photon beams propagating in 
opposite directions. Backward emission is also proposed for 
interpretation of specific features of the pulse profile (Dyks 
et al. 2005). It is suggested here that the backward propagat- 
ing waves must pass through the CFLR and undergo strong 
induced scattering. As the zeroth order approximation, one 
assumes the photon beam is confined to the magnetic merid- 
ian plane and ignoring both aberration and general relativ- 
ity, one may estimate the radial distance to the crossing 
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point (tc^c) from Vc/ cos{a + t/i*) — r,/ cos{a — ipc)- This 
gives 

where ^ 1 is assumed, and 7r/2 — a ~ 3^«/2. The polar 
angles t/ic and t/), can be expressed in terms of the relevant 
radu: Vc ~ {tc/RlcV^^ and V* ~ {r,/RLcY^'^- On sub- 
stituting them into 1191 . one finds the solution rc ~ t'*/4. 
When aberration is included, the propagation path of the 
backward emission bends toward the star in the corotating 
frame. A more rigorous treatment is outlined in Appendix. 
Since for fast rotating young pulsars or fast millisecond pul- 
sars, the critical radius rcr is close to or larger than the LC 
radius, propagation of the backward emission in the CFLR 
is not allowed. For example, as shown in Sec 4.3, the Crab 
pulsar has rcr > 0.5Rlc for n" > 1. For the backward emis- 
sion not to pass through the opaque region (r < r^r), the 
emission radius must satisfy r, > Arcr > 2Rlc- This implies 
that no backward emission can pass through the CFLR. 

For a typical pulsar, the critical radius can be smaller 
than the LC radius. Pulsars satisfy this condition (vcr < 
Rlc) must have a period 

1/5 / \ -1/5 / R \ 3/10 



P > (0.1s) 



Ann 



50 



3 X 1022 K 

2/5 



/^N-V5._^. 

\lGUzJ VlO^T^ 



(20) 



Backward emission can propagate through the region rcr < 
r < Rlc, provided that the emission is produced close to 
the LC, 4rcr < r, < Rlc- We consider PSR B0950+08 as 
an example. It has a period P = 0.253 s and a magnetic 
field B « 10** T. We estimate AQTb ~ 10^'' K and rcr ~ 
O.OURLcn"^^'^ . Backward emission can propagate through 
the CFLR if the emission region is at r, > 0.056RLC'n"'^^^ . 
Note that n" can be much larger than unity. For n" — 10, 
one has r« > 0.2Rlc- Thus even for a typical pulsar, propa- 
gation of radio emission in the CFLR is strongly constrained 
and is allowed only in the region near the LC. 



5 CONCLUSIONS 

We consider the effect of induced three-wave turbulence on 
propagation of the radio emission in a pulsar magnetosphere. 
(As remarked in Sec. 2.5, this process is sometimes referred 
to, inappropriately we believe, as 'stimulated Raman' or 'in- 
duced Raman' scattering.) The low frequency wave is either 
an electrostatic plasma wave propagating along magnetic 
field lines or an Alfven wave propagating nearly along the 
magnetic field. We discuss both cases of propagation in the 
OFLR and CFLR. In the polar (OFLR) region, the induced 
processes strongly depend on the bulk velocity of the plasma 
and on the pair multiplcity (M). Since the scattering rate 
(cf. Eq. 1131 decreases for increasing 7 and M, the induced 
three- wave scattering effect on the propagation is suppressed 
due to the relativistic bulk flow (7 ^ 1) and M 1. As- 
suming a modest bulk Lorentz factor 7 ~ lO'^ and a mul- 
tiplicity M ~ 10^ for a young pulsar like the Crab pulsar, 
and 7 ~ 10* and M ~ 10 for a fast millisecond pulsar like 
PSR 19574-21, we show that propagation of radio emission 
in the OFLR is little affected by the induced three-wave 



processes. An exception is the propagation of giant pulses. 
If the giant pulses are emitted relatively close to the star, 
they are subject to very strong induced processes that pre- 
vent them from escaping. The favored region where induced 
processes can be effective is the transition region located at 
r,/7 from the source, where the photon beam propagates 
nearly perpendicular to the magnetic field in the rest frame. 
For giant pulses not to be completely destroyed by such 
processes, their emission region needs to be close to the LC. 
Th is conclusion is cons i stent with the conclusion reached 
by iRomani fc JohnstonI i200 J) from a different argument, 
based on the alignment of giant pulses with the high energy 
emission. 



Induced three- wave interactions are generally important 
for propagation of radio emission in the CFLR. Because the 
plasma in this region is stationary in the corotating frame, 
with \r]\ = 1, the processes are much more efficient than 
in the polar region. The induced three-wave processes are 
important throughout the CFLR except for pulsars with 
a relatively long period. We discuss three specific exam- 
ples where the processes are relevant. (1) For the double 
pulsar J0737-3039B, we find that the processes discussed 
here are not efficient enough to explain the observed eclipse, 
due to the very low plasma density of the CFLR of pul- 
sar B which has a relatively long period. However, the in- 
duced three-wave scattering can be effective if one assumes 
a density much higher than the GJ density (by a factor of 
10*). It is worth noting that a much higher plasma density 
is neede d in the synchr otron /cyclotron absorption model 
as well iLvutikov fc Thompson 2005: Rafikov & GoldrcichI 
200^. (2) Our result imposes a strong constraint on the loca- 
tion of the emission region of pulsars with a pulse profile with 
widely separated leading and trailing components. There is a 
particular radial range where the trailing component may be 
swept into the CLFR due to aberration. (3) Our result limits 
the visibility of backward emission. The backward emission 
must propagate through the CFLR where it can be dispersed 
due to three wave interactions. For young or fast millisec- 
ond pulsars, the radius of the opaque region (where induced 
three-wave interactions are important) is comparable with 
or larger than the LC radius. In this case, backward emis- 
sion may not be visible. For long period pulsars, the plasma 
density in the CFLR is low and the radius of the opaque 
region is smaller than the LC. The backward emission can 
propagate through the magnetosphere provided that it is 
produced sufficiently close to the LC and that its propaga- 
tion path is outside the opaque sphere. 



An important approximation made in our discussion is 
the strong magnetic field limit. Although such an approxi- 
mation excludes the possibility of cyclotron resonance, our 
result should be valid for a finite, strong magnetic field. This 
is because the relevant frequency of the low frequency waves 
(Langmuir and Alfven modes) in the induced three-wave in- 
teractions considered here is much lower than the cyclotron 
frequency throughout the pulsar magnetosphere. 
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APPENDIX A: PHOTON BEAM IN THE 
COROTATING FRAME 

In the non-rotating frame, a photon beam is assumed to be 
emitted at in the field hne direction and to propa- 

gate along a straight line described by 



sin((T — (t) = ±r, sin(CT -I- tpt 



(All) 



x{s) = a;, ± s cos a, 



y{s) = y, ± ssina, 



(Al) 



where s = ct/Ri^c = is the distance propagated in time 
t, a is the slope of the photon beam, and ± correspond 
respectively to the forward and backward emission. Here, all 
the relevant radii, distances are in units of i?LC- We consider 
an orthogonal rotator with a dipole field given by 

X = sin^ ^, y = ("^ cosf/jsin^ (A2) 

where 

(A3) 

sm ^pB 

is a parameter that specifies a particular field line, where 
tpd is the half-opening angle of the polar cap, tpB is the 
colatitude angle. The open field lines correspond to tpB ^ 
ipd, that is C ^ 1- For a nonrotating dipole, one must have 



tan (J = 



3 cos ipt 



3 sin ip* cos Tpt ^ ^ 

To include the rotation, we consider a corotating frame 
in which the propagation path is described by 

X{s) = J?(s)sin[1'(s)], 

Y(s) = J?(s)cos[1'(s)], (A5) 

where R{s) = [X(s)2 + Y{s)^]^^^ = r(s), *(s) = i/;(s) + s, 
and 

SLrctan[x{s) /y{s)]. 



r[s 
iP{s 



(A6) 



The colatitude angle, radial distance of the emission point 
are 'I'(O) = = V(0) = tp,, R{0) = r(0) = r,. The tangent 
is 

dY 
dX 



We use 
r'(s) : 

^'( 



y'(s) _ r' cos{i) + s) - r{l + sin(V' -I- s) 
X'{s) ^ r' sin(i/) + s) + r{l + tp') cos(i/' + s) 
sin((T — s) =p r sin(i/' -I- s) 
cos((T — s) ± r cos('!/j -|- s) 

a;(s)a;'(s) +y{s)y'{s) 
r{s) 

x{s) COS cr -I- y{s) sin a 
r(s) 

3:'{s)y{s) - y'{s)x{s) 



(A7) 



= ± 



= ±sin((7-|- V), (A8) 



= ± 



r(s)2 

l/(s) COS cr — x{s) sin cr 



r(s)2 



cos(o- -I- Ip) 
r{s) 



At s = we have 



dV 
dX 



sin CT =F r, sin tpt 
cos a ikr, cos ■)/'* 



(A9) 



(AlO) 



Since the derivation of IIAIOII is in the approximation to the 
order (3 = v/c, it can also be derived by a Galilean transfor- 
mation. Assuming dY / dX\B^o = tana, i.e. the tangent of 
the beam at s = in the corotating frame, we obtain 



The -f sign corresponds to forward emission and — to back- 
ward emission. Due to the aberration effect, the photon 
beam must make an angle a ^ a to the field line in the 
observer's frame so that it is tangent to the field line in the 
corotating frame. We assume the magnetic field is a static 
dipole in the corotating frame: 



X 



C sin^ *, 



Y 



cos sin^ ^, 



(A12) 



where C, = sin 'I'd/ sin VPs. The tangent to the field line at 
s = is 



tan cr — 



3cos^*, - 1 
3 sin vp* cos 



(A13) 



In the small angle approximation <C 1, we have 7r/2 — cr ; 
3?/'*/2 and then lAlU reduces to 



cr — o" ~ ±r, . 



(A14) 



The photon beam direction at r* in the observer's frame 
advances in the rotation direction by r*. 

For the trailing component, the condition for the for- 
ward propagating photon beam to cross the closed field 
lines with ^ 1 is 7? = sin'^ ^, where R{s) = r{s) = 
[x{sf + y{sf]^^^, that is 



-I- ± 2r«s sin(7/;» -I- a] 



P''' = C'sin'(V' + s), (A15) 



[(■ 

where tantp = (r» sin ip, + s cos a)/(r* cos tp,. + s sin cr). Note 
that the — sign corresponds to the backward emission prop- 
agating into the CLFR on the leading side. The condition 
is therefore that there exists at least a nontrivial solution 
to Eq. 1A15II for s and C ^ 1 within the LC. In the case 
of forward emission, on approximating the left-hand side of 
Eq. (lAlSt by r, -I- s and using R{s) — r{s), one obtains the 
approximate form (1181 after reversing them back to their 
dimensional form r, r,/RLc, s — > s/Rlc- Similarly, the 
condition HA 1511 for the backward emission (the — sign) can 
be written approximately as r, — s ~ (i/; + s)^ in the small 
angle approximation, where t/j « (r« — 3stf},/2) / (r» — s). One 
finds that the photon beam enters the CFLR on the leading 
side at a radius rc ~ r*/4. 
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